In Silico Design of gRNA for CRISPR/Cas9-Mediated Gene Knockout.

CRISPR/Cas9 stands as a revolutionary and versatile gene editing technology. At its core, the Cas9 DNA endonuclease is guided with precision by a specifically designed single-guide RNA (gRNA). This guidance system facilitates the introduction of double-stranded breaks (DSBs) within the DNA. Subsequent imprecise repairs, mainly through the non-homologous end-joining (NHEJ) pathway, yield insertions or deletions, resulting in frameshift mutations. These mutations are instrumental in achieving the successful knockout of the target gene. In this chapter, we describe all necessary steps to create and design a gRNA for a gene knockout to a target gene before to transfer it to a target plant.

Comparing the performance of two-stage residual inclusion methods when using physician's prescribing preference as an instrumental variable: unmeasured confounding and noncollapsibility.

Aim: The first objective is to compare the performance of two-stage residual inclusion (2SRI), two-stage least square (2SLS) with the multivariable generalized linear model (GLM) in terms of the reducing unmeasured confounding bias. The second objective is to demonstrate the ability of 2SRI and 2SPS in alleviating unmeasured confounding when noncollapsibility exists. Materials & methods: This study comprises a simulation study and an empirical example from a real-world UK population health dataset (Clinical Practice Research Datalink). The instrumental variable (IV) used is based on physicians' prescribing preferences (defined by prescribing history). Results: The percent bias of 2SRI in terms of treatment effect estimates to be lower than GLM and 2SPS and was less than 15% in most scenarios. Further, 2SRI was found to be robust to mild noncollapsibility with the percent bias less than 50%. As the level of unmeasured confounding increased, the ability to alleviate the noncollapsibility decreased. Strong IVs tended to be more robust to noncollapsibility than weak IVs. Conclusion: 2SRI tends to be less biased than GLM and 2SPS in terms of estimating treatment effect. It can be robust to noncollapsibility in the case of the mild unmeasured confounding effect.

Local intraspecific aggregation in phytoplankton model communities: spatial scales of occurrence and implications for coexistence.

The coexistence of multiple phytoplankton species despite their reliance on similar resources is often explained with mean-field models assuming mixed populations. In reality, observations of phytoplankton indicate spatial aggregation at all scales, including at the scale of a few individuals. Local spatial aggregation can hinder competitive exclusion since individuals then interact mostly with other individuals of their own species, rather than competitors from different species. To evaluate how microscale spatial aggregation might explain phytoplankton diversity maintenance, an individual-based, multispecies representation of cells in a hydrodynamic environment is required. We formulate a three-dimensional and multispecies individual-based model of phytoplankton population dynamics at the Kolmogorov scale. The model is studied through both simulations and the derivation of spatial moment equations, in connection with point process theory. The spatial moment equations show a good match between theory and simulations. We parameterized the model based on phytoplankters' ecological and physical characteristics, for both large and small phytoplankton. Defining a zone of potential interactions as the overlap between nutrient depletion volumes, we show that local species composition-within the range of possible interactions-depends on the size class of phytoplankton. In small phytoplankton, individuals remain in mostly monospecific clusters. Spatial structure therefore favours intra- over inter-specific interactions for small phytoplankton, contributing to coexistence. Large phytoplankton cell neighbourhoods appear more mixed. Although some small-scale self-organizing spatial structure remains and could influence coexistence mechanisms, other factors may need to be explored to explain diversity maintenance in large phytoplankton.

Flow in temporally and spatially varying porous media: a model for transport of interstitial fluid in the brain.

Flow in a porous medium can be driven by the deformations of the boundaries of the porous domain. Such boundary deformations locally change the volume fraction accessible by the fluid, creating non-uniform porosity and permeability throughout the medium. In this work, we construct a deformation-driven porous medium transport model with spatially and temporally varying porosity and permeability that are dependent on the boundary deformations imposed on the medium. We use this model to study the transport of interstitial fluid along the basement membranes in the arterial walls of the brain. The basement membrane is modeled as a deforming annular porous channel with the compressible pore space filled with an incompressible, Newtonian fluid. The role of a forward propagating peristaltic heart pulse wave and a reverse smooth muscle contraction wave on the flow within the basement membranes is investigated. Our results identify combinations of wave amplitudes that can induce either forward or reverse transport along these transport pathways in the brain. The magnitude and direction of fluid transport predicted by our model can help in understanding the clearance of fluids and solutes along the Intramural Periarterial Drainage route and the pathology of cerebral amyloid angiopathy.

A Genuinely Hybrid, Multiscale 3D Cancer Invasion and Metastasis Modelling Framework.

We introduce in this paper substantial enhancements to a previously proposed hybrid multiscale cancer invasion modelling framework to better reflect the biological reality and dynamics of cancer. These model updates contribute to a more accurate representation of cancer dynamics, they provide deeper insights and enhance our predictive capabilities. Key updates include the integration of porous medium-like diffusion for the evolution of Epithelial-like Cancer Cells and other essential cellular constituents of the system, more realistic modelling of Epithelial-Mesenchymal Transition and Mesenchymal-Epithelial Transition models with the inclusion of Transforming Growth Factor beta within the tumour microenvironment, and the introduction of Compound Poisson Process in the Stochastic Differential Equations that describe the migration behaviour of the Mesenchymal-like Cancer Cells. Another innovative feature of the model is its extension into a multi-organ metastatic framework. This framework connects various organs through a circulatory network, enabling the study of how cancer cells spread to secondary sites.

Joint Bayesian estimation of cell dependence and gene associations in spatially resolved transcriptomic data.

Recent technologies such as spatial transcriptomics, enable the measurement of gene expressions at the single-cell level along with the spatial locations of these cells in the tissue. Spatial clustering of the cells provides valuable insights into the understanding of the functional organization of the tissue. However, most such clustering methods involve some dimension reduction that leads to a loss of the inherent dependency structure among genes at any spatial location in the tissue. This destroys valuable insights of gene co-expression patterns apart from possibly impacting spatial clustering performance. In spatial transcriptomics, the matrix-variate gene expression data, along with spatial coordinates of the single cells, provides information on both gene expression dependencies and cell spatial dependencies through its row and column covariances. In this work, we propose a joint Bayesian approach to simultaneously estimate these gene and spatial cell correlations. These estimates provide data summaries for downstream analyses. We illustrate our method with simulations and analysis of several real spatial transcriptomic datasets. Our work elucidates gene co-expression networks as well as clear spatial clustering patterns of the cells. Furthermore, our analysis reveals that downstream spatial-differential analysis may aid in the discovery of unknown cell types from known marker genes.

A comprehensive stroke risk assessment by combining atrial computational fluid dynamics simulations and functional patient data.

Stroke, a major global health concern often rooted in cardiac dynamics, demands precise risk evaluation for targeted intervention. Current risk models, like the CHA 2 DS 2 -VASc score, often lack the granularity required for personalized predictions. In this study, we present a nuanced and thorough stroke risk assessment by integrating functional insights from cardiac magnetic resonance (CMR) with patient-specific computational fluid dynamics (CFD) simulations. Our cohort, evenly split between control and stroke groups, comprises eight patients. Utilizing CINE CMR, we compute kinematic features, revealing smaller left atrial volumes for stroke patients. The incorporation of patient-specific atrial displacement into our hemodynamic simulations unveils the influence of atrial compliance on the flow fields, emphasizing the importance of LA motion in CFD simulations and challenging the conventional rigid wall assumption in hemodynamics models. Standardizing hemodynamic features with functional metrics enhances the differentiation between stroke and control cases. While standalone assessments provide limited clarity, the synergistic fusion of CMR-derived functional data and patient-informed CFD simulations offers a personalized and mechanistic understanding, distinctly segregating stroke from control cases. Specifically, our investigation reveals a crucial clinical insight: normalizing hemodynamic features based on ejection fraction fails to differentiate between stroke and control patients. Differently, when normalized with stroke volume, a clear and clinically significant distinction emerges and this holds true for both the left atrium and its appendage, providing valuable implications for precise stroke risk assessment in clinical settings. This work introduces a novel framework for seamlessly integrating hemodynamic and functional metrics, laying the groundwork for improved predictive models, and highlighting the significance of motion-informed, personalized risk assessments.

Optimizing the releasing strategy used for the biological control of the sugarcane borer Diatraea saccharalis by Trichogramma galloi with computer modeling and simulation.

One of the challenges in augmentative biological control programs is the definition of releasing strategy for natural enemies, especially when macro-organisms are involved. Important information about the density of insects to be released and frequency of releases usually requires a great number of experiments, which implies time and space that are not always readily available. In order to provide science-based responses for these questions, computational models offer an in silico option to simulate different biocontrol agent releasing scenarios. This allows decision-makers to focus their efforts to more feasible options. The major insect pest in sugarcane crops is the sugarcane borer Diatraea saccharalis, which can be managed using the egg parasitoid Trichogramma galloi. The current strategy consists in releasing 50,000 insects per hectare for each release, in three weekly releases. Here, we present a simulation model to check whether this releasing strategy is optimal against the sugarcane borer. A sensitive analysis revealed that the population of the pest is more affected by the number of releases rather than by the density of parasitoids released. Only the number of releases demonstrated an ability to drive the population curve of the pest towards a negative growth. For example, releasing a total of 600,000 insects per hectare in three releases led to a lower pest control efficacy that releasing only 250,000 insects per hectare in five releases. A higher number of releases covers a wider range of time, increasing the likelihood of releasing parasitoids at the correct time given that the egg stage is short. Based on these results, it is suggested that, if modifications to the releasing strategy are desired, increasing the number of releases from 3 to 5 at weekly intervals is most likely preferable.

In silico analysis of the wheat BBX gene family and identification of candidate genes for seed dormancy and germination.

BACKGROUND: B-box (BBX) proteins are a type of zinc finger proteins containing one or two B-box domains. They play important roles in development and diverse stress responses of plants, yet their roles in wheat remain unclear. RESULTS: In this study, 96 BBX genes were identified in the wheat genome and classified into five subfamilies. Subcellular localization prediction results showed that 68 TaBBXs were localized in the nucleus. Protein interaction prediction analysis indicated that interaction was one way that these proteins exerted their functions. Promoter analysis indicated that TaBBXs may play important roles in light signal, hormone, and stress responses. qRT-PCR analysis revealed that 14 TaBBXs were highly expressed in seeds compared with other tissues. These were probably involved in seed dormancy and germination, and their expression patterns were investigated during dormancy acquisition and release in the seeds of wheat varieties Jing 411 and Hongmangchun 21, showing significant differences in seed dormancy and germination phenotypes. Subcellular localization analysis confirmed that the three candidates TaBBX2-2 A, TaBBX4-2 A, and TaBBX11-2D were nuclear proteins. Transcriptional self-activation experiments further demonstrated that TaBBX4-2A was transcriptionally active, but TaBBX2-2A and TaBBX11-2D were not. Protein interaction analysis revealed that TaBBX2-2A, TaBBX4-2A, and TaBBX11-2D had no interaction with each other, while TaBBX2-2A and TaBBX11-2D interacted with each other, indicating that TaBBX4-2A may regulate seed dormancy and germination by transcriptional regulation, and TaBBX2-2A and TaBBX11-2D may regulate seed dormancy and germination by forming a homologous complex. CONCLUSIONS: In this study, the wheat BBX gene family was identified and characterized at the genomic level by bioinformatics analysis. These observations provide a theoretical basis for future studies on the functions of BBXs in wheat and other species.

Powerful and accurate detection of temporal gene expression patterns from multi-sample multi-stage single-cell transcriptomics data with TDEseq.

We present a non-parametric statistical method called TDEseq that takes full advantage of smoothing splines basis functions to account for the dependence of multiple time points in scRNA-seq studies, and uses hierarchical structure linear additive mixed models to model the correlated cells within an individual. As a result, TDEseq demonstrates powerful performance in identifying four potential temporal expression patterns within a specific cell type. Extensive simulation studies and the analysis of four published scRNA-seq datasets show that TDEseq can produce well-calibrated p-values and up to 20% power gain over the existing methods for detecting temporal gene expression patterns.

Collaborative virtual reality environment in disaster medicine: moving from single player to multiple learners.

BACKGROUND: The use of virtual reality (VR) in healthcare education is on the increase. In disaster medicine, it could be a solution to the cost and logistic constraints for a "full-scale" scenarios. However, VR is mainly designed for single players, which is not appropriate for the objectives pursued in disaster medicine. We decided to evaluate the educational value of using individual VR simulation in disaster medicine on a group of learners. METHODS: The VR scenario used was a reproduction of a major train crash, with 21 victims and whose objectives were START triage and first aid techniques. The sessions were carried out in multi-participant groups with different roles (active and immersed with headset, paper triage without headset, and active for communications not immersed in the headset). Their perceived self-efficacy was assessed before (T0), after (T1) and 2 months (T2) after the training. Satisfaction and confidence in learning were also measured. RESULTS: The median levels of satisfaction and confidence in learning were of 21/25 and 32/40 respectively. Their perceived self-efficacy increased significantly between T0 and T1 (p < 0.001), and remained stable until T2. The different roles of participant showed no difference in terms of satisfaction, confidence in learning or changes in perceived self-efficacy. One third of the participants agreed that the number of participants had interfered with their learning. A significant negative correlation (rS = -0.51, p = 0.002) was found between satisfaction and the fact of having been hindered by the number of participants. Around 90% of participants found the activity entertaining and found the new technologies appropriate for learning technical skills. CONCLUSIONS: This first experience of VR in a group setting is satisfactory and shows its positive effects. The limitations highlighted here will enable areas of improvement to be identified for the use of VR in disaster medicine, pending the development of multi-player tools. It would now be appropriate to analyse the impact of this type of simulation on learning and its retention over time.

spVC for the detection and interpretation of spatial gene expression variation.

Spatially resolved transcriptomics technologies have opened new avenues for understanding gene expression heterogeneity in spatial contexts. However, existing methods for identifying spatially variable genes often focus solely on statistical significance, limiting their ability to capture continuous expression patterns and integrate spot-level covariates. To address these challenges, we introduce spVC, a statistical method based on a generalized Poisson model. spVC seamlessly integrates constant and spatially varying effects of covariates, facilitating comprehensive exploration of gene expression variability and enhancing interpretability. Simulation and real data applications confirm spVC's accuracy in these tasks, highlighting its versatility in spatial transcriptomics analysis.

Future land use prediction and optimization strategy of Zhejiang Greater Bay Area coupled with ecological security multi-scenario pattern.

The land use changes driven by human activities press a incredible menace to zonal ecological security. As the most active urban cluster, the uncontrolled expansion of cities in the bay area exerts enormous pressure on the ecosystem. Therefore, from the perspective of ecological conservation, exploring future land use optimization patterns and spatial structure is extremely essential for the long-term thriving of the bay area. On this basis, this research integrated the System Dynamics model (SD) as the quantity forecast model and the PLUS model as the spatial emulation model and established the Land Use/Cover Change (LUCC) Simulation Framework by setting the constraints of Ecological Security Multi-Scenario Patterns (ESMP). By setting four scenarios in future, that is, Business As Usual (BAU), Priority of Ecological Protection (PEP), Balanced Development Scenario (BD), and Priority of Urban development (PUD), this research predicts LUCC in the Zhejiang Greater Bay Area (ZGBA) in 2035 and explored land use optimization patterns. The results indicate that by 2035, under the scenarios of BAU, BD, and PUD, the construction land will observably grow by 38.86%, 19.63%, and 83.90%, respectively, distributed mainly around the Hangzhou Bay Area, Taizhou Bay Area, and Wenzhou Bay Area, primarily achieved by sacrificing ecologically sensitive lands such as forests to achieve regional high economic growth. Under PEP, the growth of construction land retards, and forest experiences net growth (11.27%), with better landscape connectivity and more cohesive patches compared to other scenarios. Combining regional planning and analysis at the city scale, Hangzhou Bay area (Hangzhou, Huzhou, Jiaxing, Shaoxing, Ningbo) can adopt the BD development scenario, while Zhoushan, Taizhou, Wenzhou and Fuyang County of Hangzhou can adopt the PEP development scenario. This research furnishes a novel mechanism for optimizing land use pattern in ecological security perspective and offers scientific guidance for land resource management and spatial planning in ZGBA.

Adversarial AI applied to cross-user inter-domain and intra-domain adaptation in human activity recognition using wireless signals.

In recent years, researchers have successfully recognised human activities using commercially available WiFi (Wireless Fidelity) devices. The channel state information (CSI) can be gathered at the access point with the help of a network interface controller (NIC card). These CSI streams are sensitive to human body motions and produce abrupt changes (fluctuations) in their magnitude and phase values when a moving object interacts with a transmitter and receiver pair. This sensing methodology is gaining popularity compared to traditional approaches involving wearable technology, as it is a contactless sensing strategy with no cumbersome sensing equipments fitted on the target with preserved privacy since no personal information of the subject is collected. In previous investigations, internal validation statistics have been promising. However, external validation results have been poor, due to model application to varying subjects with remarkably different environments. To address this problem, we propose an adversarial Artificial Intelligence AI model that learns and utilises domain-invariant features. We analyse model results in terms of suitability for inter-domain and intra-domain alignment techniques, to identify which is better at robustly matching the source to target domain, and hence improve recognition accuracy in cross-user conditions for HAR using wireless signals. We evaluate our model performance on different target training data percentages to assess model reliability on data scarcity. After extensive evaluation, our architecture shows improved predictive performance across target training data proportions when compared to a non-adversarial model for nine cross-user conditions with comparatively less simulation time. We conclude that inter-domain alignment is preferable for HAR applications using wireless signals, and confirm that the dataset used is suitable for investigations of this type. Our architecture can form the basis of future studies using other datasets and/or investigating combined cross-environmental and cross-user features.

G Protein Activation Occurs via a Largely Universal Mechanism.

Understanding how signaling proteins like G proteins are allosterically activated is a long-standing challenge with significant biological and medical implications. Because it is difficult to directly observe such dynamic processes, much of our understanding is based on inferences from a limited number of static snapshots of relevant protein structures, mutagenesis data, and patterns of sequence conservation. Here, we use computer simulations to directly interrogate allosteric coupling in six G protein α-subunit isoforms covering all four G protein families. To analyze this data, we introduce automated methods for inferring allosteric networks from simulation data and assessing how allostery is conserved or diverged among related protein isoforms. We find that the allosteric networks in these six G protein α subunits are largely conserved and consist of two pathways, which we call pathway-I and pathway-II. This analysis predicts that pathway-I is generally dominant over pathway-II, which we experimentally corroborate by showing that mutations to pathway-I perturb nucleotide exchange more than mutations to pathway-II. In the future, insights into unique elements of each G protein family could inform the design of isoform-specific drugs. More broadly, our tools should also be useful for studying allostery in other proteins and assessing the extent to which this allostery is conserved in related proteins.

DeepDynaForecast: Phylogenetic-informed graph deep learning for epidemic transmission dynamic prediction.

In the midst of an outbreak or sustained epidemic, reliable prediction of transmission risks and patterns of spread is critical to inform public health programs. Projections of transmission growth or decline among specific risk groups can aid in optimizing interventions, particularly when resources are limited. Phylogenetic trees have been widely used in the detection of transmission chains and high-risk populations. Moreover, tree topology and the incorporation of population parameters (phylodynamics) can be useful in reconstructing the evolutionary dynamics of an epidemic across space and time among individuals. We now demonstrate the utility of phylodynamic trees for transmission modeling and forecasting, developing a phylogeny-based deep learning system, referred to as DeepDynaForecast. Our approach leverages a primal-dual graph learning structure with shortcut multi-layer aggregation, which is suited for the early identification and prediction of transmission dynamics in emerging high-risk groups. We demonstrate the accuracy of DeepDynaForecast using simulated outbreak data and the utility of the learned model using empirical, large-scale data from the human immunodeficiency virus epidemic in Florida between 2012 and 2020. Our framework is available as open-source software (MIT license) at github.com/lab-smile/DeepDynaForcast.

Perturbation Variability Does Not Influence Implicit Sensorimotor Adaptation.

Implicit adaptation has been regarded as a rigid process that automatically operates in response to movement errors to keep the sensorimotor system precisely calibrated. This hypothesis has been challenged by recent evidence suggesting flexibility in this learning process. One compelling line of evidence comes from work suggesting that this form of learning is context-dependent, with the rate of learning modulated by error history. Specifically, learning was attenuated in the presence of perturbations exhibiting high variance compared to when the perturbation is fixed. However, these findings are confounded by the fact that the adaptation system corrects for errors of different magnitudes in a non-linear manner, with the adaptive response increasing in a proportional manner to small errors and saturating to large errors. Through simulations, we show that this non-linear motor correction function is sufficient to explain the effect of perturbation variance without referring to an experience-dependent change in error sensitivity. Moreover, by controlling the distribution of errors experienced during training, we provide empirical evidence showing that there is no measurable effect of perturbation variance on implicit adaptation. As such, we argue that the evidence to date remains consistent with the rigidity assumption.

Photobleaching FRET-FLIM-ICS for quaternary structure quantification on cells. Theory and simulations.

The oligomerization of proteins is an important biological control mechanism and has several functions in activity and stability of enzymes, structural proteins, ion channels and transcription factors. The determination of the relevant oligomeric states in terms of geometry (spatial extent), oligomer size (monomer or dimer or oligomer) and affinity (amounts of monomer, dimer and oligomer) is a challenging biophysical problem. Förster resonance energy transfer and fluorescence fluctuation spectroscopy are powerful tools that are sensitive to proximity and oligomerization respectively. Here it is proposed to combine image-based lifetime-detected Forster resonance energy transfer with image correlation spectroscopy and photobleaching to determine distances, oligomer sizes and oligomer distributions. Simulations for simple oligomeric forms illustrate the potential to improve the discrimination between different quaternary states in the cellular milieu.

Forward problem of electrocardiography based on cardiac source vector orientations.

To localize the unusual cardiac activities non-invasively, one has to build a prior forward model that relates the heart, torso, and detectors. This model has to be constructed to mathematically relate the geometrical and functional activities of the heart. Several methods are available to model the prior sources in the forward problem, which results in the lead field matrix generation. In the conventional technique, the lead field assumed the fixed prior sources, and the source vector orientations were presumed to be parallel to the detector plane with the unit strength in all directions. However, the anomalies cannot always be expected to occur in the same location and orientation, leading to misinterpretation and misdiagnosis. To overcome this, the work proposes a new forward model constructed using the VCG signals of the same subject. Furthermore, three transformation methods were used to extract VCG in constructing the time-varying lead fields to steer to the orientation of the source rather than just reconstructing its activities in the inverse problem. In addition, the unit VCG loop of the acute ischemia patient was extracted to observe the changes compared to the normal subject. The abnormality condition was achieved by delaying the depolarization time by 15ms. The results involving the unit vectors of VCG demonstrated the anisotropic nature of cardiac source orientations, providing information about the heart's electrical activity.

An efficient Bayesian observer model of attractive and repulsive temporal context effects when perceiving multistable dot lattices.

In multistable dot lattices, the orientation we perceive is attracted toward the orientation we perceived in the immediately preceding stimulus and repelled from the orientation for which most evidence was present previously (Van Geert, Moors, Haaf, & Wagemans, 2022). Theoretically-inspired models have been proposed to explain the co-occurrence of attractive and repulsive context effects in multistable dot lattice tasks, but these models artificially induced an influence of the previous trial on the current one without detailing the process underlying such an influence (Gepshtein & Kubovy, 2005; Schwiedrzik et al., 2014). We conducted a simulation study to test whether the observed attractive and repulsive context effects could be explained with an efficient Bayesian observer model (Wei & Stocker, 2015). This model assumes variable encoding precision of orientations in line with their frequency of occurrence (i.e., efficient encoding) and takes the dissimilarity between stimulus space and sensory space into account. An efficient Bayesian observer model including both a stimulus and a perceptual level was needed to explain the co-occurrence of both attractive and repulsive temporal context effects. Furthermore, this model could reproduce the empirically observed strong positive correlation between individuals' attractive and repulsive effects (Van Geert et al., 2022), by assuming a positive correlation between temporal integration constants at the stimulus and the perceptual level. To conclude, the study brings evidence that efficient encoding and likelihood repulsion on the stimulus level can explain the repulsive context effect, whereas perceptual prior attraction can explain the attractive temporal context effect when perceiving multistable dot lattices.